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We study the impurity entanglement entropy Se in quantum impurity models that feature a 
Kondo-destruction quantum critical point (QCP) arising from a pseudogap in the conduction-band 
density of states or from coupling to a bosonic bath. On the local-moment (Kondo-destroyed) side 
of the QCP, the entanglement entropy contains a critical component that can be related to the 
order parameter characterizing the quantum phase transition. In Kondo models describing a spin- 
Simp, Se assumes its maximal value of ln(2Simp + 1) at the QCP and throughout the Kondo phase, 
independent of features such as particle-hole symmetry and under- or over-screening. In Anderson 
models, 5e is nonuniversal at the QCP, and at particle-hole symmetry, rises monotonically on passage 
from the local-moment phase to the Kondo phase; breaking this symmetry can lead to a cusp peak 
in Se due to a divergent charge susceptibility at the QCP. Implications of these results for quantum 
critical systems and quantum dots are discussed. 

PACS numbers: 71.10.Hf, 71.27.-|-a, 75.20.Hr 


I. INTRODUCTION 

Entanglement entropy has emerged as a powerful tool 
for quantifying correlations in pure and mixed states 
of quantum many-body systems, and particularly for 
characterizing unconventional ground states. In cases 
where a local order parameter is not sufficient, the en¬ 
tanglement entropy Se [defined in Eq. (IT|) below] can be 
used to identify nontrivial topological order [I|. It has 
been established [Ij that Se exhibits nontrivial scaling in 
the vicinity of continuous zero-temperature phase transi¬ 
tions, i.e., at quantum critical points Q. Erom an infor¬ 
mation perspective, entanglement is a generic feature of 
quantum-mechanical systems that may be used for stor¬ 
ing information as well as for quantum computation [^. 

Due to the nonlocality of entanglement entropy, there 
are relatively few interacting problems where it has 
proved possible to accurately determine Se, particularly 
in the vicinity of quantum phase transitions. In this 
context, quantum impurity problems have the advantage 
that they can be solved reliably via a number of well- 
established techniques. Conformal field theory the 
density-matrix renormalization group [H , and the numer¬ 
ical renormalization group @ have all been used to calcu¬ 
late entanglement entropy in the Kondo and spin-boson 
models. Quantum impurity models therefore provide 
a natural setting for systematic study of the interplay 
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between quantum criticality and entanglement entropy. 
Eor example, it has been established usiM the numerical 
renormalization-group (NRG) method [aT^ ^hat Se ex¬ 
hibits a cusp peak at the quantum critical point (QCP) 
that separates the localized and delocalized phases of the 
sub-ohmic spin-boson model. Whether or not a cusp 
peak in Se at the QCP is a generic feature of continu¬ 
ous impurity quantum phase transitions is an open and 
intriguing question. 

The Kondo and Anderson impurity models 0 , origi¬ 
nally shown to describe the many-body screening of the 
magnetic moments of dilute magnetic impurities in non¬ 
magnetic metals, have now found applications in other 
contexts, including but not limited to the physics of 
heavy /-electron systems and transport through quan¬ 
tum dots. Here, we study variants of these models that 
exhibit critical destruction of the Kondo effect, where 
Kondo screening is suppressed at a second-order quan¬ 
tum phase arising from the presence of a pseudogap in 
the conduction-band density of states around the Fermi 
energy [iMl and/or an additional coupling of the im¬ 
purity to a bosonic environment [l9l - [^ . Kondo destruc¬ 
tion QCPs in the Bose-Fermi Kondo model have been 
proposed to describe anomalous quantum criticality in 
heavy fermion metals [2^, and non-equilibrium crit¬ 
icality in a quantum dot with ferromagnetic leads [^ . 
Pseudogap Kondo and Anderson models have been used 
to study nonmagnetic impurities in d-wave superconduc¬ 
tors [23 and in graphene [s^. It has also been shown 
that tuning a double-quantum-dot system can produce 
a pseudogap in the effective density of states 
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While each of these models have been studied in the past, 
little is known about their entanglement entropy proper¬ 
ties, a gap that this work aims to hll. 

By combining analytic and NRG calculations, we es¬ 
tablish that the entanglement entropy between a mag¬ 
netic impurity and its environment contains a critical 
component in the vicinity of these Kondo-destruction 
QCPs. In Kondo models with a spin-S'imp impurity 
moment, we show that Se takes its maximal value of 
ln(25'imp -I- 1) at the QCP and throughout the Kondo 
phase, and decreases in a power-law fashion on entry 
into the Kondo-destroyed or local-moment phase. These 
behaviors highlight some differences between the Kondo 
and spin-boson models. In pseudogap Anderson mod¬ 
els, we show that charge fluctuations lead Se to take a 
nonuniversal value at the quantum phase transition. At 
particle-hole symmetry, Se always increases with increas¬ 
ing impurity-band hybridization P and is therefore never 
peaked at the the quantum phase transition. Away from 
particle-hole symmetry, Se can vary non-monotonically 
with P, and in some cases exhibits a cusp peak at the 
QCP. 

The remainder of this paper is organized as follows: 
Section El reviews general characteristics of the entan¬ 
glement entropy and summarizes the universal behaviors 
that we find at Kondo-destruction QCPs. Detailed anal¬ 
ysis of various Kondo models and Anderson models is 
presented in Secs. Ell and US respectively. We discuss 
our results in Sec. El and conclude in Sec. ED 


II. GENERAL CONSIDERATIONS 

Entanglement entropy captures the degree of quantum 
nonlocality in the ground-state wave function. Specifi¬ 
cally, it is a property associated with a partition of the 
system into two regions A and B that effectively “cuts” 
the ground state along the boundary between the regions. 
Upon tracing the system’s density operator p over region 
B, one obtains the reduced density operator in region A, 
Pa = Trs p. Similarly, one can trace over region A to 
obtain ps = Tr^i p. The entanglement entropy is the von 
Neumann entropy of pA or ps, i.e., 

SeiA\B) = -TvAiPAAvpA) = -TrB(/5Bhi;OB), (1) 

which measures the extent to which region A is entangled 
with region B. 

In quantum impurity problems, the entanglement en¬ 
tropy between the impurity and the rest of the system 
is defined by taking region A to contain solely the impu¬ 
rity degrees of freedom, while region B describes the host 
(i.e., the rest of the system), as shown schematically in 
Fig-E Upon tracing out the host, we obtain the impurity 
reduced density operator pimp acting in a vector space of 
dimension dimp- Equation o then gives the impurity 



FIG. 1. (Color online) Schematic representation of the model 
Hamiltonians considered in this work. The impurity degrees 
of freedom Himp are coupled via Rhost-imp (wavy line) to the 
bath degrees of freedom Rhost. We trace out region B, and 
determine the entanglement entropy of region A. 

entanglement entropy in 

<^imp 

= - ^pjnpi, (2) 

where {pi} is the set of eigenvalues of pimp. The cor¬ 
responding eigenstates {|i)} must respect the system’s 
symmetries, a constraint that allows the eigenvalues pi to 
be expressed in terms of expectation values of impurity 
operators that can readily be calculated using the NRG. 
Since the host degrees of freedom have been completely 
traced out, the impurity entanglement entropy measures 
only the entanglement between the impurity and the host 
as a whole. Details of the host—such as the number, dis¬ 
persion, and any internal interactions of the conduction 
bands and/or the bosonic baths—influence Se only in¬ 
sofar as they affect the impurity matrix elements that 
determine the eigenvalues of pimp. 

Eor a ground state of product form |^/)) = |(/))imp C 
lx)host, one can choose pi = 1 and pi = 0 for all i > 1, 
implying that Se = 0. At the other extreme, a state 
of maximal entanglement between the impurity and its 
host is described by pi = 1/dimp for all i, leading to 

e — hi dimp. 

A complication arises if the system is not in a pure 
state, as is likely to be the case when there is ground-state 
degeneracy. For example, in the trivial limit where the 
impurity and the host are decoupled, n-fold degeneracy of 
the impurity ground state results in pimp having n values 
Pi = \/n and dimp — n values pi = 0, implying that 
Se = Inn. In order to avoid such misleading indications 
of entanglement, it is necessary to break the ground-state 
degeneracy of the impurity to obtain a pure state. 

In the present work, where we treat magnetic impu¬ 
rities, the ground-state degeneracy can be lifted by the 
application of an infinitesimal local magnetic field ft-ioc 
that couples solely to the impurity through a Hamilto¬ 
nian term where Simp is the impurity spin op- 
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FIG. 2. (Color online) Order parameter Mioc(A,/iioc = 0'^) 
and spontaneous-symmetry-breaking part of the entangle¬ 
ment entropy SSe{A) vs distance A — {J — Jc)IJc from the 
quantum critical point of the Bose-Fermi Kondo model with 
an Ising-symmetry coupling Kog ~ 7.5 between an impurity 
spin Simp = I and a bosonic bath with exponent s = 0.8. On 
approach to the QCP from the local-moment phase (A < 0), 
SSe vanishes like 

erator and the Lande g factor and the Bohr magneton 
have both been set to unity. For this reason, we consider 
a two-parameter function S'e(x,/iioc), where x is a non- 
thermal, nonmagnetic parameter that tunes the system 
through a QCP at x = Xc- In many cases, we employ a 
reduced variable A = (x — Xc)/xc such that the QCP is 
located at A = 0. It also proves convenient to define the 
local-field-dependent part of the entanglement entropy 
(note the sign) 

SSe{x, hioc) = Seix, 0) - S'e(x, /lloc), (3) 

and to introduce the shorthand notations 

S+(x) = Se(x,/lloc = 0+), (4a) 

SSe{x) = SSeix,hioc = 0^) (4b) 

representing, respectively, the degeneracy-lifted entangle¬ 
ment entropy and the reduction in entanglement entropy 
due to spontaneous symmetry breaking. 

In Secs. uni and El we report results for the entangle¬ 
ment entropy in several quantum impurity Hamiltonians 
of the general form 

FI = FFhost 4” FFimp 4" FFhost -imp (5) 

where Hhost describes one or more fermionic bands as well 
as, possibly, a bosonic bath. The term iJimp describes the 
isolated impurity, and Hhost-imp accounts for the coupling 
between the host and the impurity. The fermionic bands 
are assumed to have a dispersion Ck giving rise to an 
idealized density of states 

Pc(e) = N-^ Y. = Po\e/D\^Q{D - |e|), (6) 

k 


where Nk is number of distinct k values, 0(x) is the 
Heaviside function, D is the half-bandwidth, and r is 
a band exponent. QCPs arise in pure-fermionic Kondo 
and Anderson models in cases 0 < r < 1 describing pseu- 
dogapped hosts in which pde) vanishes in a sub-linear 
fashion at the Fermi energy = 0. 

Another route to impurity QCPs is competition be¬ 
tween fermionic and bosonic environments for control of 
an impurity degree of freedom. Below we consider mod¬ 
els with a metallic fermionic density of states described 
by Eq. m with r = 0 as well as a bosonic bath assumed 
to have a dispersion Wq that gives rise to a density of 
states 

p^{uj) = Af^"^^(5(w-Wq) = A:oWc“®|w|®0(a;)0(a;c-a;), 

q 

(7) 

where Nq is the number of distinct q modes, lOc is an 
ultra-violet frequency cutoff, and s is a bath exponent 
assumed to lie in the range 1/2 < s < 1. 

In all of the models we consider it has been estab¬ 
lished that a QCP separates a Kondo phase (correspond¬ 
ing to tuning parameter values A > 0) from a local- 
moment phase (spanning A < 0), in which the Kondo 
effect is destroyed. The models include ones in which 
the Kondo phase exhibits exact screening, over-screening, 
or under-screening of the impurity spin. In each case 
an appropriate order parameter for the quantum phase 
transition is the liioc —t O'*" limit of the local magneti¬ 
zation Mioc(A,/ lloc) = ~ l™T->.o(5'4jjp), which vanishes 
throughout the Kondo phase, and in the local-moment 
phase close to the QCP obeys 

Mioc(A,/iioc = 0+)cx(-A)^, (8) 

where P is the order-parameter exponent. At the critical 
value of the tuning parameter, the local magnetization 
satisfies 

Mioc(A = 0, /lloc) OC l/llocT^'^, (9) 

where 5 is another critical exponent. 

We show very generally—independent of features such 
as particle-hole symmetry or asymmetry, the degree of 
impurity charge fluctuations, whether the Kondo phase 
involves exact, over-, or under-screening, and the pres¬ 
ence or absence of competition between fermionic bands 
and bosonic baths—that upon approach to the QCP from 
the local-moment side, the entanglement entropy satisfies 

5S,{AMoc) = ( 10 ) 

where a is a constant that depends on details of the 
model. This connection is illustrated in Fig. [51 which 
shows the variation of Mioc(A,/iioc = 0'^) and (5S'e(A) 
with reduced Kondo coupling A = (J — Jc)/Jc in the 
Ising-symmetry Bose-Fermi Kondo model [defined in 
Eqs. (fTTll below], for which we show that a = 2. 
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When combined with Eqs. (HJ and dH), Eq. OT implies 
that 

55e(A<0,hioc = 0+)oc(-A)^% (11a) 

55e(A = 0,hioc)cx Ihiocr^^', (11b) 

with 

Pe = 2/3, (12a) 

l/<5e = 2/5, (12b) 

scaling relations that are demonstrated explicitly in the 
numerical results presented below. 

We solve the quantum impurity problems introduced 
above using the NRG [H, [s^ as adapted to treat 
quantum impurity problems involving a pseudogapped 
fermionic density of states fl^ [T^ and/or a sub-ohmic 
bath of bosons jH, [H, [H, l36j| . The impurity entangle¬ 
ment entropy is found via Eq. m using reduced den¬ 
sity matrix eigenvalues pi expressed in terms of expec¬ 
tation values (converged for large NRG iteration num¬ 
bers corresponding to asymptotically low temperatures) 
of certain impurity operators specified in the sections 
that follow. We use a Wilson discretization parameter 
3 < A < 9, a range that has been shown previously 
to p rovide an accurate account of the critical exponents 
[13, [23 - E3 |. Close to the QCP, we find it necessary to 
employ quadruple-precision floating-point calculations in 
order to accurately resolve the entanglement entropy, and 
in particular, the value of 5Se{A) defined in Eq. (j4b|. 


'S'imp = 5 and channel numbers K = 1 (exactly screened 
impurity) and K = 2 (overscreened impurity). We also 
consider the (one-channel, S’imp = 5 ) Ising-symmetry 
Bose-Fermi Kondo model described by Eq. ([S]) with 

Hhost — E ^kcrcK ^kcra +E- q</q/'q> (14a) 

k,cr,a q 

iSimp = /ilocSfmp (14b) 

-flhost-imp = J Simp ’ ^ ^ 2 ^Oa'a 

<T,(7' ,a 

+ 5'5'imp ^^EK + '(’-q)’ (14c) 

(T,a V ^ q 

Here, a QCP may be present both for a metallic (r = 0) 
and a semimetallic (r > 0) conduction band. Section 
IIIIBI treats a spin-one impurity, focusing on the under¬ 
screened K = 1 pseudogap model. In each of these cases, 
we consider W (and where present g) to be held fixed and 
define the distance from criticality to be A = (J—Jc )/ Jc- 
Even with the addition of a degeneracy-lifting field 
that couples to the Kondo and Bose-Fermi Kondo 
Hamiltonians [Eqs. (TT^ and (1141) ] exhibit spin-rotation 
symmetry about the z axis and hence conserve the z 
component of total spin. This ensures that the eigen¬ 
states of Pimp can be chosen to be the conventional basis 
states |m) such that S’fj^pjm) = mjm). Then, 

s 

Se = - Pmlnpm- (15) 

m— — S 


III. KONDO AND BOSE-FERMI KONDO 
MODELS 


A. 5imp = ^ Kondo and Bose-Fermi Kondo models 


The Kondo models under consideration are described 
by Hamiltonians of the form of Eq. m with 


H\iost — ^ Ck , 

k..cr,a 

^imp — ^loc'5*impi 

3?host-imp = <4 Simp ’ ^Oaa 2 Cga'a 

,OL 

+ »'E ^Ocra ^OcrcK i 

(7,0; 


(13a) 

(13b) 


(13c) 


For a spin-1/2 Kondo impurity and in the presence of 
spin-rotation symmetry about the z axis, the eigenval¬ 
ues of Pimp are just the impurity spin-up and spin-down 
occupation probabilities 


P±l/2 = ^ ± Mioc, 

(16) 

and Eq. ([2]) reduces to 


iSe = S2(^ + Alloc), 

(17) 

where 



where Cko-a destroys a conduction electron of wave vec¬ 
tor k, spin z component a = ±1 = f,/,, channel in¬ 
dex a £ {1,... ,K}, and energy Ck satisfying Eq. (O; 
cocra = Aj, Ckcra destroys an electron of spin z 

component a and channel index a at the impurity site; 
and Simp is the spin operator for a spin-5'imp impurity. 
The host-impurity coupling is characterized by an anti¬ 
ferromagnetic exchange J > 0 and a potential scattering 
W. 

In Sec. IIII Al we consider pseudogap Kondo models 
described by Eqs. and (1131) with an impurity spin 


5 '2(a;) = —cclncc — (1 — cc) ln(l — x) (18) 

is the binary entropy function. Expanding Eq. (ITTII for 
I Alloc I < ^ gives 

Se:^ln2-2m-(4/3)M,t, (19) 

a result that holds for any S' = 1 Kondo model, irre¬ 
spective of the number and dispersion of the conduction 
bands. 

Equation (1191) implies that even in the presence of an 
infinitesimal magnetic field, the entanglement entropy 
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FIG. 3. (Color online) Entanglement entropy in the one- 
channel, S'imp = 5 pseudogap Kondo model for various com¬ 
binations of the band exponent r and the dimensionless po¬ 
tential scattering pqW. (a) Degeneracy-lifted entanglement 
entropy St vs dimensionless Kondo coupling poJ. Results 
for r = 0.2 and r = 0.3 are at particle-hole symmetry 
{W = 0). Particle-hole-asymmetric results are shown for 
r — 0.4 {poW — 0.109), r = 0.6 {poW = 0.54), and r = 0.8 
{poW = 6.2, plotting poJ/6 on the horizontal axis). In each 
case, takes its maximum value of In 2 ~ 0.693 for all 
J > Jc. (b) Spontaneous-symmetry-breaking part of the en¬ 
tanglement entropy 5S'e(A) and order-parameter Mioc (plot¬ 
ted as 2Mi^^) vs distance |A| = (Jc — J')!Jc from the QCP in 
the local-moment phase for several of the cases shown in (a). 
The coincidence of the two data sets for each r confirms Eq. 
(I20II and straight-line fits yield exponents listed in Table H] 

takes its maximum possible value of In 2 at auy mag¬ 
netic QCP and throughout the Kondo phase. This is true 
both when the impurity moment is exactly screened with 
Fermi-liquid excitations (as is the case for iC = 1) and 
when it is over-screened with a non-Fermi liquid many- 
body spectrum (as for K > 2). Taking into account also 
Eq. dHI), one finds that on approach to the QCP from the 
local-moment side (A —>• 0“), 

SSe{A) 2Ml, cx (-A)2^ (20) 

realizing Eq. OT with a = 2 as well as exemplifying Eqs. 
(Illal) and (ll2aD . 

Another important conclusion that can be drawn from 
Eqs. © and ([T^ is that 5+ vanishes only for |Mioc| —f 5 , 


corresponding to a vanishing Kondo coupling J . Even 
though the Kondo effect is destroyed for 0 < J < Jc 
and the impurity contributes In 2 to the zero-temperature 
entropy, the impurity remains entangled with its envi¬ 
ronment and the ground state cannot be represented in 
direct-product form. 

We supplement these general theoretical consider¬ 
ations with NRG results for several specific cases. 
For iSinip = 51 we only need to compute Mjoc = 
— limT^o(<5'imp)i from which Se can be obtained exactly 
via Eq. (fT71) . Figure EDa) shows the degeneracy-lifted en¬ 
tanglement entropy S'^ vs po J for the one-channel pseu¬ 
dogap Kondo model. This model is known [I^ to have 
(1) a particle-hole-symmetric Kondo-destruction QCP at 
J = Jc,s{r), W = 0 for any band exponent r on the range 
0 < r < i, and (2) a pair of particle-hole-asymmetric 
QCPs at J = Jc,a{r), W = ±Wc(r) for any r satisfy¬ 
ing r*~3/8<r<I. The figure shows, both for the 
symmetric case W = 0 (illustrated for r = 0.2 and 0.3) 
and the asymmetric one W ^ 0 (represented by data for 
r = 0.4, 0.6, and 0.8), that rises from 0 at J = 0 to 
In 2 at J = Jc, and remains pinned at In 2 throughout 
the Kondo phase. 

Figure |3|b) is a log-log plot of and 5Se vs |A| in 
the local-moment phase. Fitting to the power-law forms 
of Eqs. [5] and (Illal) yields the exponents /3 and /3e listed 
in Table m The values of /3, along with ones of 1/(5 also 
listed, improve the accuracy of those published previously 
[l3l| . For r = 0.3, 0.4, and 0.6, the scaling Eq. (I12al) 
is obeyed to within the (small) estimated nonsystematic 
error of the exponents. For r = 0.2, W = 0 and for 
r = 0.8, FF ^ 0 we find minor deviations from scaling (at 
the 0.5% and 1.5% levels, respectively). In these cases, 
the exponent 1/d takes small values that make it very 
difficult to reach values of |Mioc| <C ^ and 6Se In 2 
even at J = Jc when using the smallest magnetic fields 
that we can apply numerically (of order hioc = 1Q~^^D 
for quadruple-precision floating-point arithmetic). This 
forces estimation of (3 and /?e at J values further from Jc 
where the asymptotic power-law form is obeyed less well. 

For most values of r the entanglement exponent satis¬ 
fies ,5e < 1, meaning that dS'/' /dA diverges on approach 
to the QCP from the local-moment side. However, for 
0.35 r < i at particle-hole symmetry, as well as for 
0.35 ^ r < 0.45 away from this symmetry (exemplified 
in Table H] by r = 0.4, pqW = 0.109), one instead finds 
/3e > 1, leading to a much weaker feature in S'/" vs A at 
the location of the quantum phase transition. 

We have also performed NRG calculations for the two- 
channel pseudogap Kondo model [the case of Eq. (fT51) 
with K = 2], which is predicted [3) [il^ [13 for 0 < 
r < Cmax (where 0.23 < r^ax < 0.24) to have a particle- 
hole-symmetric QCP between a local-moment phase and 
a non-Fermi liquid overscreened-Kondo phase. Equation 
(I 20 I) is predicted to hold independent of K, and indeed 
our numerical results summarized in Table El support 
this picture: any differences between the cases K = 1 
and K = 2 arise just from the two models having their 
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r 

poW 

P 

1/5 

P. 

0.2 

0 

0.16025(1) 

0.02645(4) 

0.324(2) 

0.3 

0 

0.35499(2) 

0.07398(2) 

0.71007(4) 

0.4 

0.109 

0.57553(2) 

0.15606(3) 

1.15106(2) 

0.6 

0.54 

0.18759(2) 

0.11696(5) 

0.3754(3) 

0.8 

6.2 

0.07578(2) 

0.06373(3) 

0.156(2) 


TABLE I. Exponents /?, S, and /Je defined in Eqs. ([Hll, ([^. 
and (lllall . respectively, for the 5imp = one-channel Kondo 
model with the five combinations of the band exponent r and 
dimensionless potential scattering poW illustrated in Fig. [3] 
Parentheses enclose the estimated nonsystematic error in the 
last digit. The exponents obey Eq. (I12all within the errors, 
apart from weak violations for r = 0.2 and r = 0.8, where the 
small values of 1/5 impede accurate evaluation of P and Pe. 

r P 1/5 Pe 

0.05 0.050(1) 0.001815(2) 0.140(5) 

0.1 0.150(2) 0.013075(2) 0.327(4) 

0.15 0.330(1) 0.036233(8) 0.667(3) 

0.2 0.908(3) 0.09205(2) 1.815(9) 

TABLE 11. Exponents P, 1/5, and /3e defined in Eqs. ([5 )) . (P |l . 
and (lllall . respectively, for the S'imp = two-channel Kondo 
model at particle-symmetry [W = 0) and four different values 
of the band exponent r. Parentheses enclose the estimated 
nonsystematic error in the last digit. The exponents obey 
Eq. (|12a|) to within a margin that is within nonsystematic 
errors for r = 0.2, but becomes larger for lower values of r, 
likely due to the small values of exponent 1/5. 


own order-parameter exponents P{r). 

The final S'imp = \ model that we have studied numeri¬ 
cally is the Ising-symmetry Bose-Fermi Kondo model [Eq. 
da]- As exemplified in Fig.|2]for the case of conduction- 
band exponent r = 0 and a sub-Ohmic bosonic-bath 
exponent s = 0.8, NRG calculations once again con¬ 
firm that 8Se obeys Eq. ([2(1 on approach to the Kondo- 
destruction QCP from the local-moment side (A < 0), 
whereas in the Kondo-screened phase (A > 0) we find 
Se(A, 0) = In2. Just as in pure-fermionic Kondo mod¬ 
els, S+ rises on approach to phase boundary from the 
local-moment side to attain is maximum possible value 
at the QCP and throughout the Kondo phase. 

The behavior of the entanglement entropy in the Bose- 
Eermi Kondo model is in marked contrast with that of 
the sub-Ohmic spin-boson model, described by Eq. ([5]) 
with 


F^host — ^ ( Wq ij^q (/fq, 

n 

( 21 a) 

FFimp = ^lociS'imp ~ 

( 21 b) 

FFhost-imp = <7>5'imp + '('-q) ’ 

( 21 c) 


q 


where the dispersion Wq satisfies Eq. 0 with 0 < s < 1. 
NRG studies [ 3 , @| show that Se rises on approach to the 
QCP from either phase, and exhibits a cusp peak at a 
non-universal value that falls short of the maximum pos¬ 


sible value In 2. This departure from the entanglement 
found in the Bose-Fermi Kondo model seems surprising, 
given that both models map [H, to a classical Ising 
model with long-range ferromagnetic interactions that 
decay for large separations d as where s is the 

bosonic bath exponent. The two quantum-mechanical 
models are thought to share the same critical exponents 
at least for i < s < 1 . 

However, the classical mapping and the resulting crit¬ 
ical exponents describing the response to a local mag¬ 
netic field along the z axis take no account of the dif¬ 
ferent global symmetries of the two models. The Kondo 
model exhibits a global U{1) spin symmetry leading to 
a conserved total spin 2 component, as well as an emer¬ 
gent SU (2) spin symmetry at the Kondo-screened fixed 
point. By contrast, the spin-boson model exhibits no 
spin symmetry at the Hamiltonian level (at least in non¬ 
trivial cases where A^, and g are both nonzero), while the 
two stable fixed points exhibit emergent U{1) symmetries 
leading to conservation of total spin x at the delocalized 
fixed point (analogous the the Kondo phase) and con¬ 
servation of total spin z at the localized (local-moment) 
fixed point. 

In the absence of a globally conserved total spin z 
component we cannot assume that the eigenstates of 
the impurity reduced density operator are also eigen¬ 
states of 5';^p. Instead, the properties Trpimp = 1 and 
Tr(/5inipO’) = 2(S) are sufficient to fully specify the im¬ 
purity reduced density operator as 

Amp = 1 + (Simp) ■ cr, (22) 

where I is the 2x2 identity operator. Then the en¬ 
tanglement entropy is given exactly by Eq. ([^l with two 
eigenvalues @ 

i ± v/(^imp)^ + (^imp)^ + (^fmp)^- (23) 

Glearly, Se attains its maximum possible value of In 2 
if and only if (Simp) = 0, a condition that is satisfied 
throughout the Kondo phase of the Bose-Eermi Kondo 
model and also at the QGP. By contrast, the nonzero 
value of Ax necessary to produce a quantum phase tran¬ 
sition in the sub-Ohmic spin-boson model leads to a 
nonzero (^j^p) in both phases. We expect that for fixed 
g, Se{hioc = 0 ) decreases smoothly from In 2 as A^, in¬ 
creases from zero, and Se{hioc = 0 ) shows no feature 
at the QCP. The addition of an infinitesimal field along 
the z axis creates a nonzero (Sj^p) only in the delocal¬ 
ized phase, thereby yielding a correction SSe obeying 
Eqs. (fini) . (Illa|) . and (I12a|) . and leading 5'+ to exhibit a 
cusp peak in all cases where 2/3 < 1. We speculate that 
similar behavior would arise in the Ising-symmetry Bose- 
Fermi Kondo model in the presence of a transverse local 
field [dll , although we have not explicitly tested this con¬ 
jecture. 
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B. S'imp = 1 single-channel pseudogap Kondo model 

For a spin-1 Kondo impurity and in the presence of 
spin-rotation symmetry about the z axis, the eigenvalues 
of Pinip can be parameterized 

P±i = ^{no± Mioc), (24a) 

po = 1 - no, (24b) 

where we have introduced uq = ((‘S’fj^p)^). For such an 
impurity, we focus exclusively on the one-channel pseudo¬ 
gap Kondo model. At particle-hole symmetry {W = 0), 
this model has a QCP for any band exponent r on the 
range 0 < r < Tmax where 0.26 < r^ax < 0.27, while 
ioi W ^ 0 there is a pair of asymmetric QCPs for any 
r* ~ 0.245 < r < 1 [l2|. Each QCP separates (1) a 
local-moment phase with a decoupled spin-one impurity 
degree of freedom and a ground-state moment-squared 
= <S'imp(*S'iinp -b 1) = 2 from ( 2 ) an under-screened 
Kondo phase in which the impurity spin and the conduc¬ 
tion band form a many-body ground state with a residual 
moment-squared = 3(2 -b r )/8 for IF = 0 or = 3/4 
for IF 0. 

In the Kondo phase, the system shows an infinitesimal 
response to an infinitesimal degeneracy-lifting field, and 
we can take the eigenvalues of pimp to be = po = 1/3, 
which corresponds to M\oc = 0, no = 2/3, and 5'e(A, 0) = 
In 3. Interestingly, even though the spin-1 impurity is 
only partially screened, the entanglement entropy still 
takes its maximal value of ln( 2 S'inip + 1 ). 

In the local-moment phase, by contrast, the ground 
state exhibits spontaneously broken SU{2) symmetry, 
and it is natural that uq should rise above 2/3. How¬ 
ever, it is plausible (and we confirm below) that close to 
the QCP, deviations of no from 2/3 will be smaller than 
those of the order parameter Mioc from zero, and that the 
former can safely be neglected. More specifically, one can 
conjecture that 


no(/iioc = 0+)-2/3oc (-A)'^", (25) 

with Pn > /?. Substituting Eqs. (l24l) into Eq. setting 
no = 2/3 and keeping only leading terms in Mioc, one 
arrives at the result 

<55e(A) « cx (-A)2^ (26) 

providing a realization of Eq. CHD with a = 3/4. As was 
the case for S'imp = 4, the predicted behavior is consistent 
with Eqs. (1 11 all and (I12al) . 

Figure 0] presents NRG results for three combinations 
of the pseudogap exponent and the dimensionless po¬ 
tential scattering: (r, poW) = (0.25, 0), (0.4, 0.5), and 
(0.6, 0.8). Here, we have computed Mioc = ~(>5'fmp) 
no = ((•S'i^p)^), then used Eqs. (0]) and (1^ to find Sg. 
Figure HKa) shows that on approach to each QCP from 
the local-moment side, Mioc, — 2/3, and 5Se vanish in 
power-law fashion according to Eqs. (O, (051) . and (IllaL 
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FIG. 4. (Color online) Approach from the local-moment side 
to the quantum critical point of the one-channel, S'lmp = 1 
pseudogap Kondo model for three combinations of the pseu¬ 
dogap exponent and the dimensionless potential scattering: 
(r, poW) — (0.25, 0), (0.4, 0.5), and (0.6, 0.8). (a) Log-log 
plot of Mioc, no — 2/3, and 5Se vs |A| for r = 0.6, with the 
linear variations exemplifying Eqs. ©, dSSl), and dllal) . re¬ 
spectively. (b) Ratio SSe/Mi^^ vs |A|. The convergence of 
the ratio to 0.75 as A —>■ 0“ (albeit at a rate that decreases 
with increasing r) confirms the validity of Eq. I®- 


respectively. The exponents extracted from plots such 
as this are summarized in Table m for the three exam¬ 
ples shown in the figure, plus the particle-hole-symmetric 
QCPs for r = 0.1, 0.15, and 0.2. In each case, the ex¬ 
ponents satisfy /3„ > 3/3, confirming the conjecture that 
deviations of hq from 2/3 can be neglected in the vicin¬ 
ity of the QCP. Equation (I12al) is also satisfied with one 
minor violation (for r = 0.25) and one more significant 
one for r = 0.6 that can be attributed to the same cause 
(a small exponent l/d) as in the one-channel S'lmp = \ 
Kondo model. 

Figure IlKb) plots the ratio 5Se/Mf^^ on the same log¬ 
arithmic |A| scale as used for Figure |4l(a). For r = 
0.25, W = 0, the ratio converges rapidly to 0.75, as 
predicted in Eq. (1261) . The convergence is rather slower 
for r = 0.4, pqW = 0.5, and markedly slower for 
r = 0.6, PqW = 0.8. This trend can be explained by the 
taking into account the leading correction to Eq. (I2S1) in 















r 

poW 

P 

Pn 

P. 

0.1 

0 

0.062(5) 

0.19(1) 

0.13(1) 

0.15 

0 

0.116(1) 

0.386(5) 

0.24(1) 

0.2 

0 

0.207(1) 

0.787(4) 

0.423(4) 

0.25 

0 

0.51(1) 

2.7(1) 

1.00(5) 

0.4 

0.5 

0.08960(5) 

0.4509(3) 

0.180(1) 

0.6 

0.8 

0.04300(5) 

0.2128(2) 

0.089(1) 


TABLE III. Exponents /3, /3n, and /3e defined in Eqs. (HSJ, 
and (lllall . respectively, for the S'imp = 1, one-channel Kondo 
model six combinations of the band exponent r and the di¬ 
mensionless potential scattering poW. Parentheses enclose 
the estimated nonsystematic error in the last digit. In each 
case, Pn > 3/3, implying that variations of no are negligible 
compared to those of Mioc in the vicinity of the quantum crit¬ 
ical point. The exponents obey Eq. (fT^ within the errors, 
apart from a weak violation for r = 0.2 and a stronger one 
for r = 0.6, where the small value 1/5 = 0.02589(2) impedes 
accurate evaluation of /3 and /3e. 


situations where Pn > 2/3, namely, 

+ (27) 

and noting in Table IIIII the decrease in the value of P 
as r increases toward 1. We can extrapolate the ratio 
5Se/My^^ for A —>• 0 to obtain 0.7499(3) and 0.746(6) for 
r = 0.4 and r = 0.6, respectively, in excellent agreement 
with the predicted value a = 3/4. We therefore take 
Figure |4l(b) as providing confirmation of the relation Eq. 

(EH). 


IV. ANDERSON MODELS 

In this section we consider non-degenerate Anderson 
impurity models characterized by an impurity Hamilto¬ 
nian 

37 imp — ^d'^d 4 ” Tld_l 4 ” ^locij^dX 

where destroys an electron of energy Cd and spin z 
component cr electron at the impurity site, Uda = d\d^ 
and nd = rid^ + Udi are impurity number operators, and 
U is the on-site Coulomb interaction. Such models dif¬ 
fer from their Kondo counterparts by allowing charge 
fluctuations on the impurity site. We will show that 
such fluctuations can signihcantly modify the behavior 
of the entanglement entropy in the vicinity of the contin¬ 
uous quantum phase transition between Kondo-screened 
and local-moment (Kondo-destroyed) phases. In par¬ 
ticular, in the presence of an infinitesimal symmetry¬ 
breaking held, the entanglement entropy in general nei¬ 
ther exhibits a peak at the QCP (as it does in the sub- 
Ohmic spin-boson model) nor attains its maximum pos¬ 
sible value at the QCP and throughout the Kondo phase 
(as is the case in the Kondo and Bose-Fermi Kondo mod¬ 
els discussed in Sec. mil- 


Under conditions of spin-rotation symmetry about the 
2 ; axis, the eigenstates of the impurity reduced density 
operator can be taken [4^ to be the conventional basis 
states | 0 ), \a) = and | 2 ) = d|d|| 0 ), with eigenval¬ 

ues that can be written as 


P 2 = («dt«d 4 .), (29a) 

Pt = //2 + Mioc, (29b) 

p^ = //2-Mioc, (29c) 

Po = l-/-P 2 , (29d) 


in terms of the local-moment fraction (i.e., single¬ 
occupation probability) / = (n^) — 2 p 2 and the local 
moment M\oc = {\{nd-\ — ''^dj,)), such that 0 < / < 1 and 
|-^ioc| < //2. Then Eq. (E]) can be written 


Se = fS2 


Ml 


loc 


+ (1 - f)S2 


P2 

1 -/ 


(30) 


which can be interpreted as the sum of a binary spin en¬ 
tanglement entropy with weight / and a binary charge 
entanglement entropy with weight (1 — /). The parallel 
between the spin and charge parts of Se can be made 
clearer by defining a “local charge” Qioc = P 2 — Po, sat¬ 
isfying IQiocI < 1 - /, so that 


P2 1 Q\qc 

w = 2 + 2(1-/)- 


(31) 


In the absence of a symmetry-breaking held hioc, we 
can set Mioc = 0 in Eqs. (l29l) . and use Eq. dJ]) to obtain 


S'e(^loc = 0) = - (1 - / - P 2 ) ln(l - f -P 2 ) 

-/ln(//2)-P 2 lnp 2 - (32) 


Since / and are nontrivial functions of U, e^, and P, it 
is clear that Eq. (15^ can encompass much richer behavior 
than the corresponding Kondo-model result 5e(hioc = 
0) = ln2. Differentiating Eq. (15^ with respect to P 2 for 
a hxed local-moment fraction / yields 


\Op 2 Jf P2 


(33) 


which implies that for a given value of /, 5'e(/iioc = 0 ) 
is greatest for equal occupation of the impurity conhgu- 
rations | 0 ) and | 2 ) (i.e., Qioc = 0 ), and is smallest when 
one or other of the conhgurations is ruled out. 

In the local-moment phase near the boundary with the 
Kondo phase, we expect an inhnitesimal held hioc = O’*" 
to establish a local magnetization |Mioc| -C 5 with neg¬ 
ligible shift of / and p 2 - Under these circumstances, an 
expansion of Se in powers of Mioc shows the spontaneous- 
symmetry-breaking part of the entanglement entropy to 
be 


55e(A):^2MiL//cx(-A)2/3, (34) 

providing yet another realization of Eqs. (ITUl) . (Illal) . and 
(I12a|) . this time with a = 2/f. Since the local-moment 
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fraction satisfies 0 < / < 1 , the critical part of the entan¬ 
glement in nondegenerate Anderson models is generally 
enhanced by a factor of 1 // compared to its counterpart 
in ■S'imp = \ Kondo models [see Eq. d^ ]. 

For the purposes of numerical study, we focus on the 
one-channel Anderson impurity model described by the 
Hamiltonian in Eq. ([5]) with Hhost as given in Eq. (I13a|) . 
i/imp as in Eq. (1^ and 

^^host-imp = H.C.) , (35) 

The hybridization matrix element V between the impu¬ 
rity site and the conduction band is conventionally re¬ 
expressed in terms of the hybridization width E = . 

We have computed {n^ and {rid i)i then used Eqs. 
@ and ( 1 ^ to find 5'e. 

In the following we take U and td to be fixed, either 
at particle-hole symmetry (e^ = —17/2) or away from it 
{td ^ ~U/2). We then find the location of the QCP at 
a critical hybridization width Tc{U, Cd), and thereafter 
define the distance from criticality as A = (E — Ec)/Ec. 
The critical responses to a local magnetic field near the 
symmetric and asymmetric QCPs of the pseudogap An¬ 
derson model belong in the same universality classes as 
the respective QCPs of the S'imp = 5 Kondo model [H- 

[ 3 , [3113. 

In the subsections that follow, we first consider two 
special cases {U = —2td and U = 00 ) in which Se{hioc = 
0) in Eq. (1321) reduces to a function of one variable, the 
local-moment fraction, thereby simplifying analysis of the 
behavior of the entanglement entropy in the vicinity of 
the QCP. Afterward, we present illustrative examples of 
the entanglement properties for more general cases. 


A. Particle-hole symmetry: U = —2ed 

For the particle-hole-symmetric case = —17/2, we 
have (rid) = 1 and pQ = p 2 = {1 — /)/2. As a result, Eq. 
(EH) reduces to 

Se{hioc = 0) = S2{f)+\n2, ( 36 ) 

which increases inonotonically from In 2 to In 4 as j / — ^ j 
decreases from i to 0. Equation (155)) tells us that for a 
given value of /, 5'e(hioc = 0 ) is greater for this symmetric 
case than for any value Cd ^ —U12 that leads to po ^ P 2 - 
In the conventional situation of on-site Coulomb re¬ 
pulsion (i.e., U >0), there is a monotonic evolution 
of the local-moment fraction from / = 1 at E = 0 to 
/ —>■ i for E 17, reflecting the increased admixture 
of the Ud = 0, 2 excited configurations into the = 1 
ground state of the isolated impurity. This behavior, 
which is exemplified in Fig.jSJa) for r = 0.4 and two cases, 
U/D = 0.005 and 0.5, leads to a inonotonic increase in 
Se{hioc = 0) from In2 at E = 0 to In4 for E ^ 17. 
Both / and Se{h\oc = 0) increase in a smooth, featureless 





FIG. 5. (Color online) Particle-hole-symmetric pseudogap 
Anderson model with band exponent r — 0.4 for different val¬ 
ues of U/D: (a) Local-moment fraction / and (b) symmetry- 
lifted entanglement entropy , both vs hybridization width 
r scaled by its critical value Ec. In the limit of large E, 
approaches its maximum value of In 4 « 1.386. Inset 
to (b): Spontaneous-symmetry-breaking part of the entan¬ 
glement entropy SSe vs |A| [where A = (E — rc)/rc) for 
the two cases shown in (a) and (b). The linear variations 
of the points on this log-log plot are consistent with Eq. (fna 
with /3e = 1.8288(1). (c) Local-field-dependent part of the 
entanglement entropy SSe vs hioc at the critical hybridiza¬ 
tion width, showing behavior consistent with Eq. (fllbl) with 
l/Se = 0.3703(1). 


fashion as E rises through Ec. [The mapping U ^ —U, 
td —>■ —id, E —> E takes one to an Anderson model with 
on-site attraction, which for a pseudogap density of states 
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features quantum phase transitions between local-charge 
and charge-Kondo phases [i^ ■ Since the mapping trans¬ 
forms / —>• 1 — /, it preserves Se{hioc = 0). We will not 
consider cases U < 0 any further in this paper.] 

We have already argued that the effect of a local mag¬ 
netic held hioc = 0 ^ is to reduce the entanglement en¬ 
tropy by an amount that varies near the QCP according 
to Eq. (IM)) . Since the QCP always occurs at a local mo¬ 
ment fraction / > ^, <S'e(/iioc = 0 ) has a positive slope at 
r = Pc- It is therefore the case that the degeneracy-lifted 
entanglement entropy has a positive slope on both sides 
of the quantum phase transition. In other words, it is 
impossible for 5+ to exhibit a peak at the QCP. In cases 
where the order-parameter exponent satishes /3 < I, we 
expect a divergence of dS^ /dP on approach to the QCP 
from the local-moment side. By contrast, for 2/3 > 1, 
the spontaneous-symmetry-breaking part of the entan¬ 
glement entropy should vanish faster than the the linear 
variation Se{hioc = 0 ) and should therefore be essen¬ 
tially featureless on passing through the quantum phase 
transition. 

Figure Eljb) plots over a broad range 0 < P/Pc < 2 
for r = 0.4, a case [3, 01 where /3 = 0.91440(2) > i. 
The curves for different values of U are quantitatively dif¬ 
ferent, but share the same principal features: a smooth 
rise of from zero in the decoupled-impurity limit 
r /U —>■ 0, with a linear variation through P = Pc [not 
shown in detail in Figure [5][b)] leading to a saturation 
5+ —>■ ln4 in the uncorrelated limit T/U —>■ oo. 

We find that for a given value of A, in the particle- 
hole-symmetric Anderson model generically exceeds 
in the counterpart S'imp = ^ Kondo model with the same 
band exponent r. This is a natural consequence of the 
Anderson impurity being entangled with both the spin 
and charge degrees of freedom of its environment. In the 
limit U/D oo, however, the form of 5'+ (A) for the An¬ 

derson model converges to the corresponding function for 
the Kondo model. [The Kondo function 5'+(A) for sev¬ 
eral values of r can be extracted from Fig.[31[a) by rescal¬ 
ing the horizontal axis from p^J to A = (J — Jc)/Jc] 
Conversely, as U/D approaches zero, the critical hy¬ 
bridization vanishes [i^ as Tc/D oc [U/DY~'' while S'+ 
becomes ever closer to In4 throughout the region T > Pc; 
a trend that smoothly merges into the physics of a non¬ 
interacting resonant level {U = Cd = 0) where 5+ = ln4 
for any P > 0 . 

Examination of the spontaneous-symmetry-breaking 
part of the entanglement entropy allows quantitative 
testing of Eq. (l34l) . The inset to Fig. [5][b) illustrates 
log-log plots of d5'e(A) that for U/D = 0.5 and 0.005 
can be fitted to Eq. (Illal) with jdg = 1.8288(1), matching 
within estimated error the scaling prediction of Eq. (I12al) 
using the aforementioned value /3 = 0.91440(2). Table HVl 
shows that the ratio is captured to high accu¬ 

racy by Eq. (IMl) . 

We conclude our discussion of the particle-hole sym¬ 
metric Anderson model by presenting results for the ef¬ 
fect of a finite local magnetic field. Figure [5Kc) plots the 


U/D /c _ 2/f, _ b 

0.005 0.637220 3.13863 3.1386(1) 

0.05 0.646976 3.09131 3.0913(2) 

0.5 0.689800 2.89939 2.8994(1) 

TABLE IV. Properties at the quantum critical point of the 
particle-hole-symmetric pseudogap Anderson model for r = 
0.4 and three values of U/D: local-moment fraction /c (from 
NRG), and the predicted and actual values of SSe/My^^ in a 
field hioc = 0^, i.e., 2//c based on Eq. dSl) and the ratio b 
computed directly from NRG values of SSe and Mioc. 

local-held-dependent part of the entanglement entropy 
i5S'e(A, hioc) as a function of hioc for fixed A = 0 (i.e., 
at the critical hybridization width). The linear variation 
of the data for three different values of U/D fulfills the 
prediction of Eq. (jllbp . while the fitting of quadruple- 
precision results spanning fields down to h\oc = 10~^^D 
(data not shown) yields an exponent 1/Se = 0.37032(2) 
in excellent agreement with the value 0.37032(4) de¬ 
duced from Eq. (|12bD using the previously known value 
1/d = 0.18516(2). 

B. Maximal particle-hole asymmetry: U = oo 

For the case U = oo oi maximal particle-hole asymme¬ 
try, the parameters entering Eq. (l32ll reduce to p 2 = 0, 
/ = (rid), and po = 1- f, so that 

5e(/iioc = 0)=^2(/) + /ln2. (37) 

For a given /, this value is smaller by (1 — /) In 2 than its 
counterpart for the symmetric model; indeed, Eq. (13311 
indicates that for fixed /, Se{h\oc = 0 ) takes its smallest 
value when U = oo (and also when U = —oo, leading to 
Po = 0). Under variation of the local-moment fraction, 
'S'e(hioc = 0) increases from 0 at / = 0 to a maximum 
value of In 3 at / = 2/3, and then decreases to reach In 2 
at / = 1 . 

For given values of r and e^, a quantum phase transi¬ 
tion between local-moment and Kondo phases will occur 
at some F = Fc. Unlike the situation at particle-hole 
symmetry, where Se{hioc = 0) varies linearly with A near 
the QCP, the variation for any ^ —U/2 is described 
by 

5e(A, 0) - 5e(0,0) ~ A|A|1-^ sgn A, (38) 

where 7 is the charge-susceptibility exponent at the QCP 
and A may be positive or negative. For U = 00 , this 
variation can be deduced by writing 

d5e(A, 0)/dA = [d5e(/*loc = 0)/d/] • (d//dA) (39) 

and noting that while dSe{hioc = 0)/df is regular near the 
QCP, df /d/s. = d{nd)/dA oc |A|“^ can be nonanalytic 
[ITI . 0. For 0.55 < r < 1, the charge susceptibility 
exponent 7 is positive [4^ , signaling a critical divergence 
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FIG. 6. (Color online) U = oo pseudogap Anderson model 
with band exponent r = 0.6 for —td/D = 0.05 and 0.5: (a) 
Local-moment fraction / and (b) degeneracy-lifted entangle¬ 
ment entropy , both vs hybridization width F scaled by its 
critical value F^. For ed = —0.511, / decreases monotonically 
with increasing F, and passes through 2/3 at F/Fc ~ 1.07, 
where St peaks at In 3. For ed = —0.05D, / drops through 
2/3 (and St peaks) at F/Fc ~ 1.00085, and / (and St) 
reaches a minimum at F /Vc — 2.05 (Inset). On the wide scale 
of F /Vc shown in this figure, the location of the peak in St is 
virtually indistinguishable from the position of the quantum 
critical point. 


of impurity charge fluctuations on approach to the QCP 
from either phase. The asymmetric QCPs for r* ~ 3/8 < 
r < 0.55 instead show a non-divergent charge response, 
i.e., 7 = 0 , a behavior also displayed at the symmetric 
QCPs that exist for = —C//2 and 0 < r < i. 

As a result of the combination of the local-moment 
contribution from Eq. (1341) and the charge contribution 
from Eq. (|38)) . the degeneracy-lifted entanglement en¬ 
tropy will exhibit a cusp peak at the QCP provided that 
( 1 ) dSe{h\ac = 0 )/(iA < 0 at A = 0 to ensure that 
dS'^/d/s. < 0 for A = 0+, and (2) 2/3 < 1 — 7 to ensure 
that dSf/dA > 0 for A = 0“. Previous NRG calcu¬ 
lations have shown that condition ( 2 ) is satisfied at the 
particle-hole-asymmetric QCPs for 0.42 < r < 1 [4^. In 
what follows, therefore, we focus on whether condition 
( 1 ) is satisfied. 

The sign of dSe{h\oc = 0)/dA at A = 0 can be de¬ 
termined using Eq. dSS- From Eq. dSZl) we see that 


dSe{hioc = 0)/df is positive for / < 2/3 and negative 
for /> 2/3. In the regime < 0 that admits interest¬ 
ing many-body physics, the local-moment fraction has 
limits / —?► 1 for P —>• 0 and / —;► i for P —)► 00 (i^ . 
It might therefore appear plausible that df /dT < 0 and 
hence df /dA < 0 for all intermediate values of P. How¬ 
ever, NRG calculations show this assumption to be cor¬ 
rect only for large values of \ed\/D. For smaller \€d\/D, 
f instead has a minimum at a finite value of P, beyond 
which it increases to approach i from below. 

We present data here solely for the representative case 
r = 0 . 6 , but have obtained qualitatively similar results 
for other r values. Figure ISDa) plots / vs P/Pc for 
ed/D = —0.05 and —0.5. For ed = —0.5D, f decreases 
monotonically with increasing hybridization width, pass¬ 
ing through 2/3 at P ~ l.OTPc where 5'e(/iioc = 0) 
rises to a smooth peak at its maximal value In 3. For 
ed = —0.05D, f drops through 2/3 and Se{hioc = 0) 
peaks at r ~ l.OOOSSPc, barely into the Kondo phase. 
However, in contrast to its behavior for ed = —0.513, 
the local-moment fraction then reaches a minimum value 
/ ~ 0.834 ~ 0.759 In 3 at P ~ 2.05rc before rising back 
toward / = i. Upon further decrease of \ed\/D (not 
shown), the peak in Se{h\oc = 0) moves ever closer to 
r = Pc and the minimum value of / and a related mini¬ 
mum in Se{hioc) = 0 become deeper while the location of 
these minima remains significantly above Pc. However, 
for all Cd < 0 it appears that dSe{h\oc = 0)/df > 0 and 
df /dA > 0 at A = 0, meaning that the conditions are 
never met for the occurrence of a peak in Sf precisely at 
the QGP. 


C. General case 

Finally we turn to cases of intermediate particle-hole 
symmetry, for which no simplification of Eq. (1321) is possi¬ 
ble. We focus once again on the case r = 0.6 representa¬ 
tive of the range of band exponents in which 2/3 < 1 — 7, 
potentially allowing for the occurrence of a peak in Sf 
at the QCP. 

Figure [7] shows results over a wide range of P/Pc val¬ 
ues for ed = —0.05Z3 and for U/D = 0.051, 0.055, and 
0.075. As U gets closer to —ed, particle-hole asym¬ 
metry grows stronger, as do changes near P = P^ in 
the local-moment fraction [Fig. EDa)], the double occu¬ 
pancy [Fig. EDb)], and the degeneracy-lifted entangle¬ 
ment entropy [Fig. [71[c)]. For U/D = 0.075 and 0.055, 
Se{hioc = 0 ) peaks just inside the local-moment phase 
at A ~ —1.2 X 10“^ and —3.2 x 10“"*^, respectively. The 
negative slope of Se{hioc = 0) at P = Pc, combined with 
the order-parameter exponent satisfying 2/3 < 1 — 7 is 
sufficient to create a cusp peak in Sf at the QCP. By 
contrast, for U/D = 0.051, S'e(hioc = 0) peaks just inside 
the Kondo phase at A ~ 5 x 10“®, thereby preventing the 
occurrence of any peak in Sf at A = 0 [49|. However, 
the tiny displacement of the maximum in Se{hioc = 0 ) 
from A = 0 means that in practice it will prove very 
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FIG. 7. (Color online) Pseudogap Anderson model with band 
exponent r = 0.6, impurity level energy = —0.05D, and 
three values oi U ^ —2ed- (a) Local-moment fraction /, (b) 
double occupancy P 2 , and (c) degeneracy-lifted entanglement 
entropy St vs hybridization width F scaled by its critical 
value Fc. A peak in St occurs precisely at the the quantum 
critical point for U/D = 0.075 and 0.055, but shifted slightly 
into the Kondo phase (F > Fc) for U = 0.051II. These cases 
are indistinguishable on the wide scale of F/Fc shown in this 
figure. 


hard to distinguish cases where there really is a peak in 
the degeneracy-lifted entanglement entropy at the QCP 
from ones where a peak lies close by inside the Kondo 
phase. 

Deeper into the Kondo phase, 5+ = S'e(/iioc = 0) ex¬ 
hibits a broad minimum, which is centered at F/Fc ~ 
1.65, 1.52, and 1.13 for U/D = 0.051, 0.055, and 0.075, 


respectively. The minimum in 5'+ weakens as U increases 
toward —2ed, consistent with the monotonic rise in en¬ 
tanglement entropy throughout the Kondo phase that is 
seen at particle-hole symmetry [Fig. [SKc)]. 

The physics near the QCP is shown in more detail in 
Fig. [51 once again for a representative case r = 0.6 and 
Cd = —0.0511. Figure ISKa) plots the critical part of the 
entanglement entropy in the absence of a local magnetic 
held for U = 0.055Z1. The divergent charge susceptibil¬ 
ity produces a nonanalytic variation of iS'e(A, 0 ) —5'e(0, 0 ) 
that is well-captured by Eq. (l38l) , even though the power- 
law variation in the local-moment phase is conhned to a 
rather narrow region of |A| values due to the aforemen¬ 
tioned peak in Se{hioc = 0) at A = —3.2 x 10“^. The 
inferred exponent 7 = 0.23(2) is fully consistent with the 
directly-computed [lEllll charge-susceptibility exponent 
7 = 0.210(2). Figure |5][b) shows the local-magnetic- 
held response at the QCP (A = 0) for U/D = 0.051, 
0.055, and 0.075. Fitting to Eq. (Illbl) yields an expo¬ 
nent l/i5e = 0.2340(1), in excellent agreement with Eq. 
(I12b|) given that 2/5 = 0.23392(8) based on Table [H 


V. DISCUSSION 

One universal feature of our results is the presence of 
a nonzero entanglement on entry to the local-moment 
(Kondo-destroyed) phase. Such a residual entanglement 
implies that the ground state is not a simple product 
of an impurity state and an environmental state. This 
result has significant implications for theoretical and nu¬ 
merical descriptions of the Kondo-destroyed phase. For 
example, within a large-A mean held theory of the pseu¬ 
dogap Kondo model 0, the local moment is repre¬ 
sented with fermionic spinons fa and the effective Hamil¬ 
tonian is a resonant-level model with a hybridization 
b = {b)MF = {fa^oa)MF (where 5 is a bosonic opera¬ 
tor). At this level, Kondo destruction corresponds to 
6 —>■ 0, implying that the local moment is completely 
free and no longer entangled with the conduction band. 
Thus, such a static mean-held theory cannot reproduce 
the nonzero entanglement entropy that we hnd in the 
Kondo-destroyed phase. Our results can be understood, 
however, in terms of a bosonic operator 5(w) that has a 
vanishing static component and give rises to a dynamical 
Kondo effect. Our results also imply that the Kondo- 
destroyed phase cannot be captured in variational quan¬ 
tum Monte Carlo studies of the Kondo lattice that treat 
6 as a static variational parameter. It will be interesting 
to try and consider more general variational wave func¬ 
tions that can treat the Kondo-destroyed phase more ac¬ 
curately. 

Kondo-destroyed quantum critical points have been in¬ 
voked to understand the unconventional quantum criti¬ 
cality observed in experiments on heavy-fermion metals 
[soil . As a result of the failure of the Hertz-Millis-Moriya 
theory of the spin-density-wave transition to de¬ 

scribe the experimental data [^. the concept of local 
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FIG. 8. (Color online) Critical behavior of the entanglement 
entropy Se for the pseudogap Anderson model with r — 0.6, 
ed = —0.05D, and U ^ —2ed'- (a) Shift in the zero-field en¬ 
tanglement entropy for U = 0.0551? on moving away from the 
quantum critical point into the local-moment phase (A < 0) 
and the Kondo phase (A). Power-law fits to Eq. (I38II yield 
7 = 0.23(1). (b) Local-field-dependent part of the entangle¬ 
ment entropy 5Se vs at the critical hybridization width, 
for U/D = 0.051, 0.055, and 0.075. The data are consistent 
with Eq. (fTTbl) with l/(5e = 0.2340(1). 


quantum criticality [2^ has been used to understand the 
energy-over-temperature scaling in the dynamic spin sus¬ 
ceptibility, the presence of an additional energy scale, 
and a jump in the Fermi-surface volume. The theory 
of local quantum criticality is based on the extended- 
dynamical mean-field theory of the Kondo lattice (^ . 
which finds that for sufficiently strong quantum fluctua¬ 
tions the Kondo effect is indeed destroyed at the antifer¬ 
romagnetic QCP. The results of the present study imply 
that a continuous loss of entanglement is expected at 
the Kondo-breakdown QCPs believed to occur in certain 
heavy-fermion systems [55l |. 

Our results might be tested using engineered realiza¬ 
tions of some of the models, including single-electron 
transistors attached to ferromagnetic leads and cer¬ 
tain double-quantum-dot systems [3l|, . Equation (15^ 

shows that in the Kondo phase, knowledge of (n^) and 
{nd^Udx) is sufficient to determine Sg- Since the charge 


on a quantum dot can be measured using a variety of 
quantum point-contact techniques [5^ , and it is now be¬ 
coming possible to probe double occupancy using opti¬ 
cal spectroscopy techniques there is a very exciting 
possibility of being able to extract the impurity entan¬ 
glement entropy directly from experiments. 

The remarkable simplicity of the variation of the im¬ 
purity entanglement entropy near a Kondo-destruction 
QCP, as embodied in Eq. m, can be attributed to the 
application of Eq. © with subsystem A containing just 
the impurity degrees of freedom and subsystem B en¬ 
compassing just the host degrees of freedom. If one con¬ 
sidered a different partition of the system in which host 
degrees of freedom were split between subsystems A and 
B, then Sg would probe entanglement within the host, 
which should be much more sensitive than impurity-host 
entanglement to details such as the number of conduction 
bands and/or the presence of bosonic baths. The critical 
behavior near impurity QCPs of entanglement entropy 
defined using alternative system partitions forms an in¬ 
teresting open question for future work. 


VI. CONCLUSIONS 

We have studied the quantum mechanical entangle¬ 
ment between a magnetic impurity and its environment 
in several models that feature critical destruction of the 
Kondo effect. In the Kondo-destroyed phase of each 
model studied, we have identified a term in the entangle¬ 
ment entropy varying with a critical exponent (5g = 2/3, 
where /3 is the critical exponent governing the order pa¬ 
rameter characterizing the quantum phase transition. In 
addition, we have established that the response of Sg 
to a local magnetic field gives rise to a part of the en¬ 
tanglement entropy that varies with a critical exponent 
\/5g = 2/(5, where 5 is the critical exponent governing 
the response of the order parameter at the QCP to a lo¬ 
cal magnetic field. We have established very generally 
that in Kondo models, the ratio of the critical part of 
the entanglement entropy to the square of the order pa¬ 
rameter depends only on the magnitude of the impurity 
spin, and not on the number of conduction channels or 
the presence of an additional bosonic bath. In nondegen¬ 
erate Anderson models, this ratio is enhanced over its 
value in the 5'imp = i Kondo model by a factor of the 
inverse of the impurity’s local-moment fraction. 

Our investigation has shown the absence of any univer¬ 
sal behavior on approach to a Kondo-destruction quan¬ 
tum critical point from the Kondo (disordered) phase. In 
all variants of the Kondo model that we have considered, 
Sg remains pinned at its maximal value of ln(25'imp + 1) 
throughout the Kondo phase. However, in Anderson 
models the presence of charge fluctuations introduces 
terms that depend on the values of {nd-\ + nd\) and 
{nd'^nd])- We have shown that in the pseudogap An¬ 
derson model with a band exponent r on the range 
0.55 r < 1, charge fluctuations produce a nonanalytic 
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leading variation of Se near the QCP with a critical ex¬ 
ponent that depends only on r. Away from particle-hole 
symmetry, Se may rise on approach to the QCP from 
the Kondo side, producing a cusp peak in Se precisely 
at the quantum phase transition. However, we also find 
situations in which the entanglement entropy decreases 
continuously, albeit nonanalytically, on passing from the 
Kondo phase to the local-moment phase. 
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